No Proissart Bound from Gluon Saturation 
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In the previous work Q , we showed that while the nonlinear QCD evolution equation of Balitsky 
and Kovchegov (BK) leads to saturation of the scattering amplitude locally in impact parameter 
space, it does not unitarize the total cross section. This result was recently challenged by Ferreiro, 
Iancu, Itakura and McLerran (FIIM) who claim that the dipole-hadron cross section computed from 
the BK equation saturates the Froissart bound. In this comment, we point to a fundamental error 
in the argument of FIIM which invalidates their conclusion. We show that if the total cross section 
violates unitarity for a coloured scattering probe, it does so also for a colourless scattering probe. 



In Ref. H, we demonstrated in two independent ways 
that the dipole-hadron cross section evolved by the BK 
equation |jj does not unitarize. The first derivation was 
given in the target rest frame where the evolution resides 
in the projectile wave function. This derivation uses ex- 
plicitly colourless scattering probes. The second deriva- 
tion was given within a formulation where the evolution 
is ascribed to the target. This part of the argument can 
be construed as calculating the scattering probability of 
a coloured probe. This latter point seems to have caused 
some confusion. In particular, the authors of Ref. 
claimed recently that our conclusion holds for coloured 
probes only and that "for the physically interesting case 
where the external probe is a (colourless) dipole, there 
is no problem with unitarity at all." The main purpose 
of this note is to show that once the Froissart bound is 
violated for a coloured probe, it is also violated for a 
colourless scattering probe. We do this by considering 
the BK evolution directly as the evolution of the target. 

In the large N c limit, the BK equation is a closed equa- 
tion for the t — In 1/x rapidity dependence of the scat- 
tering probability N(x,y) = ^Tr(l - W(x)U{y)) of a 
colour singlet dipole with charges at points x and y, 



[N(x,z) 



cPz (x — y) 2 
2n (x — z) 2 (y — z) 2 
N{y, z) - N{x, y) - N(x, z)N{z, y)} 



(1) 



where a„ = 



The unitary matrix U(x) depends on 



the colour fields in the target and represents the eikonal 
factor picked up by the fundamental charge at transverse 
coordinate x while propagating through the target field. 

In the following we concentrate on the evolution of the 
scattering probability of a dipole of fixed size r = x — y 
at fixed impact parameter b = ^p'. We will sometimes 
write N(r,b) instead of N(x,y). Schematically the x,y 
dependence of N at any given rapidity is depicted on 
Fig.l. The circle of radius Ro(t) bounds what one could 
call "the interior of the hadron". For impact parameters 
b = > Ro (t) the scattering probability is vanishingly 
small for dipoles of any size. The radius R(t, Q 2 ) bounds 
the region inside which N(r, b) oc 1 for fixed r 2 = Q~ 2 
and outside which N(r, &) < 1. More precise definitions 



of Ro(t) and R(t,Q 2 ) will be given shortly. The total 
cross section for the scattering of a dipole of size Q^ 1 on 
the target is related to R(t, Q 2 ) as 



a(t,Q 2 ) 



d 2 bN(t, 6,< 



2irR 2 (t,Q 2 



(2) 




Black region = black disk 
for a dipole of size 1/Q 

Grey region = black disk 
for a coloured probe 

White regioh^Tegion^n 
which N(x,y)<l for dipoles 
of arbitrary size. Saturation 
scale is not defined here. 

FIG. 1: The three impact parameter regions of a hadronic 
target, relevant for the total dipole-hadron cross section 

We follow the expansion of Ro(t) and R(t, Q 2 ) starting 
from some initial value of rapidity to. In order to do 
this we place the initial dipole at the impact parameter 
b > R (t ). We then follow the evolution of N(r, b). We 
show that: 

a) . Initially, as long as b > Ro(t) the main contribution 
to the growth of N(r, b) in eq.([l|) comes from the points in 
the central region \z — x\ = z — y| ~ 6 with z < Ro(t) on 
the right hand side of eq. (|l ) . This long range contribu- 
tion leads to exponential growth R^(t) oc exp{a s et}. It 
also generates a power like tail in the scattering probabil- 
ity N(r,b) oc r 2 /b 2 and a nonvanishing saturation scale 
Q s (b) in the faraway region of impact parameters. 

b) . After Ro(t) reaches the impact parameter b, fur- 
ther evolution in the vicinity of b is dominated by local 
contributions from \z — y\ ~ \z — x\ *C b. These do 
not affect the growth of Ro(t) but instead lead to a fast 
growth of Q s (b). As a result it takes only a short time 
for Q s (t, b) to reach the value Q, so that R(t, Q 2 ) grows 
exponentially in t. 
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We explain that the analysis of Ref. fails on the fol- 
lowing point. The initial dipole in Ref.Q is placed in 
the "grey" region at R(t ,Q 2 ) < b < R (t Q ). This set 
up by construction disregards the evolution in the re- 
gions of transverse plane where the growth is dominated 
by the long range contributions. It therefore ab initio 
misses the most important contributions to the growth 
of the total cross section. The crucial issue is not how 
fast the "grey" region becomes "black", but instead how 
fast the "white" region becomes "grey". By placing the 
initial dipole in the "grey" region, the authors of Ref.Q 
prevent themselves from addressing this question. 

We start by discussing some general properties of the 
scattering probability N(x,y). 

Properties of N(x, y) 

The dipole scattering probability N(x,y) depends on 
the size r = x — y of the dipole {x, y} and on its impact 
parameter b = | (x + y) with respect to the hadronic tar- 
get. For sufficiently small b and r, it may be characterized 
in terms of the saturation scale Q s (t, b) 

N(x,y) = l-Tr(l-U^x)U(y)) 

~ l-exp[-Q 2 s {t,b)(x-y) 2 ] . (3) 

It is important to keep in mind that the last equality 
only holds for small impact parameters. This last equal- 
ity implies that for sufficiently large dipole size r, the 
scattering probability always approaches unity. However 
in reality at large enough impact parameters the scatter- 
ing probability is perturbatively small for dipoles of ar- 
bitrary size. This follows since the colour fields far from 
the hadron vanish and thus U(x) ~ 1 + O(g). Hence 
N(x,y) ~ 0(a s ) <C 1. Thus, for sufficiently large impact 
parameter, there is no transverse scale 1/Q s (t,b) above 
which the dipole scattering probability is of order unity, 
in contrast to the last line in (|J). 

To clarify the limited applicability of the concept of 
saturation scale Q s {t, b) in characterizing the dipole scat- 
tering probability, we illustrate the r- and 6-dependence 
of N(x,y) in more detail: 

Consider first a small dipole r <C Ro(t) where Ro(t) 
denotes the radius within which the target colour fields 
are large at rapidity t. In this case, the two legs x and y of 
the dipole see the target field of approximately the same 
normalized strength S(b), and the scattering probability 
can thus be written in the factorized form 

N(x,y)=n(r/L(b)) S(b). (4) 

Here, L(b) characterizes the scale of correlations of the 
target field (correlation length) as a function of impact 
parameter. For impact parameter b < Ro(t) where 



target colour fields are large, the colour field correla- 
tion length is equivalent to the inverse saturation scale, 
L(b) = 1/Q s (t,b). For larger impact parameters the 
correlation scale L(b) is still a meaningful concept but 
Q s (t, b) does not exist. 

At initial rapidity to, a physically sensible choice of the 
target field strength at impact parameter b is 

S(b) = 1 for b<Ro(t ), (5) 

S(b) ~ exp[-M(6-E (io))] for b>R (t ),(6) 

where M is some mass scale. One should keep in mind, 
however, that even though this initial profile shows an ex- 
ponential fall-off at large distances, there is no guarantee 
that this functional shape is preserved in the evolution 
via (|l|). In fact, as seen in cq. (g) below, it is not. 

As for the properties of the function n in (jij) , one ex- 
pects that it is a monotonically increasing function of its 
argument. Since N(x, y) is always smaller than unity and 
since its maximal value is set by the strength of the target 
colour field, we have n(r/L(b)) < 1. For small dipoles, 
the scattering probability is vanishingly small and grows 
as r 2 . Thus even the central region of the target is trans- 
parent for sufficiently small dipoles. For larger dipoles of 
size r oc L(b) the central region is black. However, for 
the profile (||g) the last line of (||) holds only as long as 
S(b) = 1, or b < Rq. In short, the initial scattering prob- 
ability N(to, x, y) can be parameterized in terms of some 
saturation scale Q s (to,b) only for impact parameters ly- 
ing in the "grey" region of the target, at b < Ro(to). The 
same of course applies also to scattering probability at 
any given value of t. 

The situation is very different for dipoles of large size 
where the factorization (^) does not hold. This is the 
case e.g. for dipoles with r = b sitting at large impact 
parameter b ^> Ro(t). If r and b are orthogonal, then 
both legs of the dipole are in the "white" region (i.e., 
x,y 3> Ro{t)) and the scattering probability is small, 
N(x, y) <C 1. If r and b are parallel instead, then one of 
the legs is in the "grey" region (i.e., x < Ro(t)) while the 
other is in the "white" region (i.e., y ^> Ro{t), U{y) ~ 
1). As a consequence, the dipole scattering probability 
loses in this latter case any dependence on the size of the 
dipole and its impact parameter. It can be viewed as 
the scattering probability 1 — jj-Tr (U"^ (x)) of a coloured 
probe in the sense 

N(x,y) = ±-Tr(l -tf(x)U(y)) 

~ 1 - i-Tr([/t( x) ) , £ 0r y » Ro{t) . ( 7 ) 

For x < Ro(t) the target fields are large, and < U(x) >= 
0. Thus N(x,y) = 1. Accordingly, the radius Ro(t) char- 
acterizes the size of the black disk for a coloured probe. 
Again, for dipoles for which (^) applies, the scattering 
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probability N(x,y) cannot be characterized in terms of 
a saturation scale Q s . 

To sum up, the saturation scale Q s {t, b) is a meaningful 
concept only in the region of impact parameter space 
which is black for a coloured probe, b < Ro(t). For larger 
impact parameter, the scattering probability is smaller 
than unity for arbitrarily large dipoles and the last line 
of (H) does not apply. 



Long-range Contributions violate the Froissart 
Bound 



We consider now a small-size dipole r = x—y = l/Q <C 
Ro(to) which has at initial rapidity to a negligible in- 
teraction probability with the hadronic target since it 
sits at large impact parameter b 3> Ro(to)- We ask how 
rapidly its scattering probability grows to order unity as 
a function of rapidity t. We first establish that the initial 
growth of N(x,y) is dominated by long-range contribu- 
tions to the z-integral in (|l|). 

If the newly generated dipole leg z in eq.([l]) sits at 
transverse position in the neighbourhood of x, y, i.e. in 
the white region, then N(x, z) ~ N(y, z) ~ N(x, y) <C 1 
and the non-linear term in the BK equation can be ne- 
glected. As a consequence, the short-range contribu- 
tion grows proportional to the local scattering probability 
which was taken to have an exponential fall-off at initial 
rapidity, 



jN{x,y) 



at t—to 



oc a s e 



-M(b-R„(t„)) 



(8) 



short 



For the long-range contribution, a lower bound of the 
growth of N(x,y) can be obtained by considering con- 
tributions from the "grey" region only. With the initial 
dipole in the white region (i.e. N(x,y) ~ 0) but the leg 
z in the "grey" region (i.e. N(x,z) ~ 1, N(y,z) ~ 1), 
the BK equation (pi) gives 



±N{ X ,y) 



long 



(9) 



Comparison of eqs. (||) and (^) shows that at suffi- 
ciently large impact parameter b 3> Ro(t), the dominant 
growth of the dipole scattering probability comes from 
long-range contributions. 

It is important that R Q (t) appearing on the right hand 
side of @ is the radius of the black region for a coloured 
probe in the precise sense specified in (Q) . Namely, while 
the right hand side of ([!]) depends on colourless dipoles 
{x,y}, {x,z} and {y, z}, only one of their coloured legs 
z touches the target. Thus, the scattering probability 
for these dipoles depends only on the strength of the 
target-induced colour rotation 1 — jj-Tr {U' (z)) for this 



one coloured leg, i.e. it is the scattering probability of a 
coloured probe. 

As shown in Ref. [lj || and as rederived in Ref. [0, 
the black disk radius Ro(t) for a coloured probe grows 
exponentially with rapidity, 



R (t) = Ra(t ) exp a s e(t - t ) 



(10) 



Substituting this into eq. (0) and integrating up to rapid- 
ity 



-1m 



ae Ro(t ) ' 



we find 



N(r, b) 



r 
b 2 



(11) 



(12) 



Thus, as advertized we find that by the time the radius 
of the " grey" region Rq (t) reaches the impact parameter 
b 3> Ro(to), the scattering probability N(r,b) is not ex- 
ponentially suppressed as suggested by eqs. (^,||) but 
rather has a long power tail. 

Since we expect R(t,Q 2 ) < Ro(t), we can not use the 
approximation eq.(||) beyond the rapidity t\. We can 
nevertheless establish the exponential growth of R(t, Q 2 ). 
Starting at rapidity t\, the point b enters the "grey" re- 
gion. From this point on the evolution is dominated by 
"local" contributions. The concept of saturation scale is 
well defined at b for t > t\. As shown in many previous 
works []5| , and rederived in the saturation momentum 
grows exponentially with rapidity 



Q s (t, b) = Q.(ti, b) exp [a a A (t - ti)] 
Thus at rapidity t% such that 



1 , Q 

<2 - h = — In 7- , , 
a\ Qs(ti,o) 



(13) 



(14) 



the scattering probability N(r, b) reaches unity. 

Putting eqs.(pd|) and ( Eh together we find that the 
radius of the black disk R(t, Q 2 ) reaches a given impact 
parameter b at rapidity t 



t 



ae Ro(to) 



Q 



Qs(h,b) 



(15) 



We have not calculated the exact value of Q s (ti,b), but 
even without doing so we can put an upper bound on t. 
Whatever the value of Q s (ti,b), it is greater than the 
inverse of the impact parameter Thus substitut- 

ing b for Q~ 1 (t\, b) in eq.jlil) and adding cq.(|ll|) we find 
the upper bound on the time that it takes the black disk 
R(t, Q 2 ) to reach a given impact parameter b 



1 R(t,Q 2 ) 
ae Ro(t ) 



-^-ln[R(t,Q 2 )Q] . (16) 
a\ 
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Thus 

a(t, Q 2 ) ~ 2nR{t, Q 2 ) > 2ir[R*{t a )Q^]^ . 

(17) 

We stress that eq. ( |l7| ) is the lower bound on R(t, Q 2 ). For 
example in the case when Q s {ti) is a fixed fe-independent 
scale, the rates of growth of R(t,Q 2 ) and Ro(t) are the 



Critique of FIIM 

The mistake in the argument of Ref. j| can now be 
easily understood. The set-up of Ref. j| is in principle in- 
adequate for studying the growth of R(t, Q 2 ). In Ref. 
the saturation scale at impact parameter b is assumed to 
be defined and to be perturbatively large at initial ra- 
pidity to, Q s (to,b) 3> Aqcd- Thus, by definition, this 
set-up ignores evolution in the region b > Ro(to) where 
the target field strength is initially small and where the 
saturation scale does not exist at initial rapidity to- As is 
clear from the discussion of eq. (^]) above, the power tails 
generated in this region quickly start dominating the to- 
tal cross section, and therefore cannot be ignored. It is 
this powerlike fall off of N(b) (or equivalently S(b)) and 
not the initial exponential one, that has to be taken as 
the "basic distribution" which is further evolved by local 
effects. This has been ignored in Ref. 0. 

The physics is quite simple. We know from Refs. 
jl], pi U that Ro(t) grows exponentially with t. Then 
eq.(|lj) tells us that for R (t ) < b < Ro(t), that is up to 
distances exponentially larger than the initial radius Rq, 
the long range contributions generate a powerlike tail for 
the scattering probability. This is to be contrasted with 
the exponential fall off assumed for the initial profile. We 
can now use the same argument as in Ref. [Q to establish 
that the total cross section of a dipole of size r = 1/Q 
grows exponentially. Starting from rapidity t\, where t\ 
is given in eq. dlT|) , the evolution of N(r, b) at the impact 
parameter b is governed by the BFKL equation. The ini- 
tial condition for this BFKL evolution at t\ is furnished 
by the previous evolution from to to t\ which is dom- 
inated by the long range contributions to eq.(0). This 
initial condition is given by eq . ( |l2|) . Thus all we have to 
do is to replace in eqs.(4.5) of Rcf.|4j] the exponential tail 
expj__2m7r&] by the power law tail r 2 /b 2 = (Q 2 & 2 ) -1 (see 
eq.(|lj)) and take into account the fact that the BFKL 
evolution goes on not for the whole rapidity interval t, 
but rather only for t — t\. Neglecting the so called diffu- 
sion term as in Ref. 141 we obtain: 



N(Q 2 ,b) = exp 



\n[Q 2 b 2 ]+Lua s (t~ti)-\^jt 



R(t, Q 2 ), we obtain 
R(t,Q 2 ) = 



-] e i? (io)" 



(19) 



(18) 



Equating the scattering probability to one at b 



Here A is the saturation momentum in the central re- 
gion of the target at the initial value of rapidity. Given 
that 2 A < ui, eq.(|l^) is consistent with the lower bound 
eq. ([IT]) . Thus we see again that the radius of the black 
disk for the scattering of a colourless dipole grows ex- 
ponentially with t. The result of Refs.|l], || stands con- 
firmed. 



Final comments 

We finally give an intuitive physics argument of why 
the total cross section of a colourless dipole violates the 
Froissart bound if it is violated by a coloured scatter- 
ing probe. To this end, we first recall an argument 
due to Heisenberg: consider a theory with a mass gap, 
where the profile of the distribution of matter density 
in any target must decay exponentially at the periphery, 
p(b) oc exp [— mb]. As this target is struck by a projectile, 
in order to produce an inelastic scattering event at least 
one particle must be produced. Assuming that the scat- 
tering is local in the impact parameter plane, the region 
of the overlap of the probe and the target must therefore 
contain energy at least equal to the mass of the lightest 
particle, m. For scattering at energy E — s/m in the 
frame where all the energy resides in the target, the tar- 
get energy density is Ep(b). Thus the scattering can only 
take place for impact parameters smaller than those that 
satisfy Eexp[— mb] = m. Thus b max = — In [s/m 2 ], 
which is equivalent to the Froissart bound. Conversely, 
if the cross section grows as a power of energy, then the 
density distribution in the target is not exponential but 
power like. With p{b) oc b~ x one obtains b max oc s ~ and 
this violates the Froissart bound. Assume now that the 
total cross section for a coloured probe (a quark, say) vio- 
lates unitarity, i.e., the scattering probability for a quark 
at impact parameter b is proportional to p{b) oc y^. A 
quark-antiquark dipole of transverse size r will compare 
the target field strength at two points separated by the 
distance r. Its scattering probability is thus proportional 
to r-^p(b) oc Tp^rr, which is still a power-law fall-off. 
Hence, the total dipole-hadron cross section violates the 
Froissart bound if the total quark-hadron cross section 
does so. 



Note added 

After the appearance of the present paper, the authors 
of H added a note to [Q agreeing with our conclusions. 
Unfortunately they also attribute to the present paper as 
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well as to our previous papers Q,J|] several statements 
which we never made nor intended to make. For exam- 
ple, they attribute to our works "the statement that one 
cannot compute the high-energy cross section using the 
ideas of (perturbative) saturation" . To the contrary, it 
is the main point of Rcf. Q to explain how the ideas of 
perturbative saturation can be used for the calculation of 
high-energy cross sections. For lack of space, we refrain 
here from giving further examples of the same kind. Al- 
though Ref. Q is quite misleading in the way it presents 
some aspects of our work, we are confident that the care- 
ful reader will easily identify these inaccuracies. 
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